Constraints on nuclear matter parameters of an Effective Chiral Model 
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Within an effective non-linear chiral model, we evaluate nuclear matter parameters exploiting the 
uncertainties in the nuclear saturation properties. The model is sternly constrained with minimal 
free parameters, which display the interlink between nuclear incompressibility {K), the nucleon 
effective mass {m*), the pion decay constant (/^) and the ct— meson mass {rria)- The best fit among 
the various parameter set is then extracted and employed to study the resulting Equation of state 
(EOS). Further, we also discuss the consequences of imposing constraints on nuclear EOS from 
Heavy-Ion collision and other phenomenological model predictions. 
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I. INTRODUCTION 

The framework of Quantum Hadrodynamics [J, 0] as 
an elegant and consistent theoretical treatment of finite 
nuclei as well as infinite nuclear matter laid down the pil- 
lars of relativistic theories which seem to provide solution 
to the so called "t/ie Coester band^ problem [1, How- 
ever, our present knowledge of nuclear matter is confined 
around nuclear saturation density (po ~ 3 x lO^^gcm^'^) 
and therefore, in order to have some meaningful correla- 
tions while extrapolating to higher densities, the nuclear 
equation of state (EOS) must satisfy certain minimum 
criteria quantified as the "nuclear saturation properties" , 
which are the physical constants of nature. Basically it 
is understood that the inherited uncertainty at po gets 
more pronounced at higher densities (3 — 10 po), rele- 
vant to astrophysical context such as the modeling of 
neutron stars. In this context, the two most important 
quantities which play vital role and are known to have 
substantial impact on the EOS are the nucleon effec- 
tive mass and the nuclear incompressibility Q . Ironi- 
cally, these two properties arc not very well determined 
and they posses large uncertainty. The nuclear incom- 
pressibility derived from nuclear measurements and as- 
trophysical observations exhibit a broad range of values 
K = (180 - 800) MeV 0. Further the non-relativistic 
and the relativistic models fails to agree to a commom 
consensus. The non-relativistic calculations predict the 
coinpression moduli in the range K = (210 — 2AQ)MeV 
[1jS[i3|j whereas, relativistic calculations predicts it in 
the range (200 - 300) MeV [H [ll|. Apart from that 
we are inevitably marred by the uncertainty in the de- 
termination of mass of the scalar meson (ct- meson). The 
attractive force resulting from the scalar sector is respon- 
sible for the intermediate range attraction which, along 
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with the repulsive vector forces provides the saturation 
mechanism for nuclear matter 0] . The estimate from the 
Particle Data Group quotes the mass of this scalar me- 
son 7o(600)' or cr-meson in the range (400-1200) MeV 
p^ . A recent estimate however, for sigma meson mass 
is found to be 513 ± 32 MeV [11]. 

Phenomenologically, parallel to the well known a — uj 
model, preferably known as the Walecka model p], 1^. llH . 
chiral models [11, [13, [H, H S [HI, [13 have been devel- 
oped and were applied to nuclear matter studies. Chiral 
symmetry is a symmetry of strong interactions in the 
limit of vanishing quark masses and is desirable in any 
relativistic theory. However, because the current quark 
masses are small but finite this symmetry can be con- 
sidered as an approximate symmetry. This symmetry is 
spontaneously broken in the ground state. In the context 
of (T— models, the a— field (which carries the quantum 
numbers of the vacuum) attains a finite vacuum expec- 
tation value (cr) = ctq = /tt . Equivalently, the potential 
for the (7— field attains a minimum at /^r [23l . [2^ . The 
value of reflects the strength of the symmetry break- 
ing and experimentally it is found to be /tt ~ 131 MeV 

m 

Time and again, the aforesaid facts and figures em- 
phasize the need to address the importance of imposing 
constraints to the EOS to narrow down the uncertain- 
ties both experimentally and theoretically. Arguably, to 
address these issues, one needs a model that has the de- 
sired attributes of the relativistic framework and which 
can be successfully applied to various nuclear force prob- 
lem both in the vicinity of po a-s well as at higher densities 
with the same set of parameters. With this motivation, 
we choose a model [25| which embodies chiral symme- 
try and has minimum number of free parameters (total 
five) to reproduce the saturation properties. The sponta- 
neous breaking of chiral symmetry relates the mass of the 
hadrons to the vacuum expectation value of the scalar 
field and thus naturally restricts the parameters of the 
model. Therefore, the present study, apart from test- 
ing the reliability of the model, puts valuable constraint 
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on the EOS based on the pion decay constant and brings 
out correlations between between the pion decay constant 
(/tt); the (T— meson mass (mo-), the nuclear incompress- 
ibility (K) and the nucleon effective mass (m*). 

In section 2, we briefly describe basic ingredients of 
the hadronic model and the energy and pressure of many 
baryonic system is computed following the mean-field 
ansatz. Subsequent section (Section 3) describes the 
methodology to evaluate the model parameters. In Sec- 
tion 4, we extract the best fit among the various parame- 
ter of the model and apply it to study the resulting EOS 
of symmetric nuclear matter. In the result and discussion 
section, we discuss the consequences of imposing various 
constraints on the model parameters and finally, we con- 
clude with some important findings of this work. 

II. THE EFFECTIVE CHIRAL MODEL 

Using the cliiral sigma model with dynamically gen- 
erated mass for vector meson, Glendening studied finite 
temperature aspects of nuclear matter and its applica- 
tion to neutron stars [l^. However, there the p— meson 
and its isospin symmetry influence was not considered. 
Although a nice framework respecting chiral symmetry, 
a drawback was its unacceptable high incompressibility 
and in the subsequent extension of the model [13], the 
mass of the vector meson is not generated dynamically. 
The model that we consider j25j in our present analy- 
sis embodies higher orders of the scalar field in addi- 



tion to the dynamically generated mass of the vector me- 
son. Without higher order in scalar field interactions, the 
model was first employed to study high density matter 
[2^. To bring down the resulting high incompressibil- 
ity, non-linear interaction in the scalar field was included 
in later work [2^ and subsequently applied to study nu- 
clear matter at finite temperature [30|. The success of 
the model then motivated us to generalize it to include 
the octet of baryons and to st udy hyperon rich matter 
and properties of neutron star [25l. l3ll|. However, in ear- 
lier works, the parameter sets that were employed were 
not studied and analyzed in detail with respect to the 
inherent vacuum properties of chiral symmetry. More- 
over, rather than a phenomenological fit the parameters 
must be constrained meaningfully, so that the resulting 
EOS is more realistic and purposeful. Motivated by this, 
we presently try to explore the consequences of imposing 
stringent constraint on the model parameters not only 
with properties known at saturation density but also on 
the resulting EOS with other phenomenological model 
predictions and experimental data at high density. In ad- 
dition to that, the correlation between various quantities 
with the vacuum value of the scalar field naturally spells 
out definite interlink between them. We now proceed to 
describe the salient features of the present model. The 
effective Lagrangian of the model interacting through the 
exchange of the pseudo-scalar meson tt, the scalar meson 
a, the vector meson lo and the iso-vector p— meson is 
given by: 
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The first line of the above Lagrangian represents the 
interaction of the nucleon isospin doublet ipB with the 
aforesaid mesons. In the second line we have the kinetic 
and the non-linear terms in the pseudo-scalar-isovector 
pion field '7?', the scalar field 'cr', and higher order terms 
of the scalar field in terms of the invariant combination of 
the two i.e., x"^ = tt^+o--^. Finally in the last line, we have 
the field strength and the mass term for the vector field 
'a;' and the iso-vector field '/t? meson. ga,guj and gp are 
the usual meson-nucleon coupling strength of the scalar, 
vector and the iso-vector fields respectively. Here we shall 
be concerned only with the normal non-pion condensed 
state of matter, so we take < 7? >= and also TOtt = 0. 

The interaction of the scalar and the pseudoscalar 
mesons with the vector boson generates a dynamical mass 



for the vector bosons through spontaneous breaking of 
the chiral symmetry with scalar field attaining the vac- 
uum expectation value xq. Then the mass of the nucleon 
(m), the scalar {m.„) and the vector meson mass {m^), 
are related to xq through 

m^gaXo, rrio = V2XX0, = g^xo ■ (2) 

To obtain the equation of state, we revert to the mean- 
field procedure in which, one assumes the mesonic fields 
to be uniform i.e., without any quantum fluctuations. We 
recall here that this approach has been extensively used 
to obtain field-theoretical EoS for high density matter 
[27| . and gets increasingly valid when the source terms 
are large The details of the present model and its 
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attributes such as the derivation of the equation of mo- 
tion of the meson fields and its equation of state {e & P) 
can be found in our preceding work [H, [3l| . For the sake 
of completeness however, we write down the meson field 
equations in the mean-field ansatz. The vector field (cu), 
the scalar field (a) (in terms of Y ~ x/xo = m* /m) and 



the isovcctor field (p) is respectively given by 



Wo 
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3 PB- 
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The quantity pb and ps are the vector and the scalar 
density defined as, 



Pb 



7 



(27r)3 



(6) 



In the above, 'fcj?' is the fermi momenta of the baryon and 
7 = 4 (symmetric matter) is the spin degeneracy factor. 
For symmetric nuclear matter {N ~ Z), we neglect the 
contribution from the p— meson. The nucleon effective 
mass is then m* = Ym and = /m^ are c^j = /iriti 
are the scalar and vector coupling parameters that enters 
in our calculations. 



PS 
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The total energy density 'e' and pressure 'P' of sym- 
(7) metric nuclear matter for a given baryon density is: 
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(8) 
(9) 



The meson field equations for u (eqn 3) and cr-mcson 
(eqn. 4) are solved self-consistently at a fixed baryon 
density to obtain the respective field strengths and the 
corresponding energy density and pressure is calculated. 



III. EVALUATION OF MODEL PARAMETERS 

Having calculated the the thermodynamic quantities 
such as the energy density and the pressure, our primary 
aim is to evaluate the set of parameters for the EoS that 
satisfies the nuclear matter properties defined at normal 
nuclear matter density (po) at zero temperature. As dis- 
cussed earlier, a desirable and valid EoS must satisfy the 
saturation properties of symmetric nuclear matter and 
the parameters of the model can be adjusted to fit those. 
Similar procedure has been adopted in Ref. [H, HI] to 
evaluate the parameters of the mean-field models. 

What we have in our hand is the set of five saturation 



properties of nuclear matter that a EOS has to satisfy, 
they are the Binding energy per nucleon (sa — 16.3)MeV, 
the saturation density (po ~ 0.153 fm~^), the nuclear 
incompressibility (167 — 380)MeV, the nucleon effective 
mass (m* jm = 0.75 — 0.90) and the asymmetry energy 
coefficient ( J = 32 ± 4)McV, all defined at po, the nu- 
clear saturation density. However, it can be seen that 
the uncertainty in their values enables us to extract and 
study the parameters within the specified range or with 
the variation thereof, in order to analyze their effect on a 
particular EOS. The five parameters of the present model 
that are to be evaluated are the three meson-nucleon cou- 
pling constants (Ccr, Ct^, Cp) and the two higher order 
scalar field constants (6 & c). 

The individual contributions to the energy density for 
symmetric nuclear matter (eqn. (8)) can be abbreviated 



as, 



(10) 



4 



where, 



Ik = — ^ / k'^dky/k'^ + m*2 
27r 



(11) 



Consequently, the respective energy contributions can 
be expressed in terms of these specified values at the 
saturation density. Using eqn. (14) and eqn. (15), they 
are given as. 



Srr = 



m2(l-y2)2 



2^,3 



IGm^c^, Co- 



(1-r^) 



2n4 



-(1-n 



(12) 



and 



po(?Ti - ai) - ^(Sefc + m* ps) 



(16) 



and 



2y2 ^ 



(13) 



where pB ~ Pn + Pp is the total baryon density which 
is the sum of the neutron density 'p„' and the proton 
density 'pp'. The relative neutron excess is then given 
by 6 = [pn ~ Pp)/ pb- At the standard state pB = Po, 
the nuclear matter saturation density and ^ = 0. Con- 
sequently, the standard state is then specified by the ar- 
gument {po,0), and the energy per particle is e{po,0) = 
e / po - m = ai = -16.3 MeV for symmetric nuclear matter. 
The nuclear matter EOS derived earlier can be expressed 
in terms of the nuclear energy density e as, 



/9o(m - ai) - 2(4^*: - m*ps) 



(17) 



where Ps is the scalar density defined in eqn. (7), an- 
alytically which is given by. 



1 



Ps 



kpEp — In 
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E 



F \ *2 
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In the above equations, m* = Ym is the effective nu- 
cleon mass and Ep = y^kp + m*^ is the effective energy 
of the nuclcon carying momenta kp. 

From eqn. (13), the vector coupling {C^j) can be read- 
ily evaluated using the relation 



£k + £cr + Slj = Po{m - ai). 



(14) 



From the the equilibrium condition P((Oo,0) = 0, we 
have, 



-e + Pb 



de 
dpB 



3^ PS 



- £n 



0. (15) 



2^2 



(19) 



with given by eqn. (17), for a specified value oiY = 
m* /m defined at po. 

Similarly, using the equation of motion for the scalar 
field (eqn. (4)), the scalar coupling can be calculated 
using the relation 



mY 



(i-y2) 



-{l-Y^f + 



-{l-Y^Y + 



2c. 



^pI 



2y4 



(20) 



In the above expression, the higher order scalar field 
couplings constants 'b' and 'c' are unknown, but they 
can be solved simultaneously to obtain the respective pa- 
rameters. To compute the constants of the higher order 
scalar field, we use the equation of motion of the scalar 
field (eqn. (4)) and eqn. (12). From eqn. (12), we get 



y4(i- 
2C(T Cu 



mps 



r(i-y2)2 



c(i - r2 



(i-y2) 



(22) 



c(l-y') ^ 16£.c^c^ 2m2c^ 46 

m2c^ (1 - y2)3 (1 - + 3 ■ ^ ) 

From the equation of the motion of scalar field, we get 



Substituting eqn. (21) in eqn. (22) leads us to the ex- 
pression to calculate the higher order scalar field constant 
'6'. which is. 
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IV. RESULTS AND DISCUSSIONS 



(23) 



Similarly, the higher order constant 'c' in the scalar 
field can be computed from the relation, 



8c, 



y(i - r2)3 
(i-r2)2- 



y4(l-r2)3 (l_y2)4 



(24) 



The calculation of C^^ is straight forward, but eqn. 
(20), (23) and (24) can be solved simultaneously numeri- 
cally for a given initial values of Co-, b and c, the solution 
of which would thus return the set of values for a desired 
value of F at po- 

Finally, for studying asymmetric matter, we need to 
incorporate the effect of iso- vector p— meson and the cou- 
pling for the p— meson has to be obtained by fixing the 
asymmetry energy coefficient J w 32 ±4 MeV [U at po- 
Accordingly, the p— meson coupling constant {Cp) can 
be fixed using the relation, 



Cp k'p 
127r2 



(25) 



where Cp = gl/m'^ and kp ~ {Gtt^pb/^Y^"^- 

Thus the model parameters are evaluated solving equa- 
tions (19), (20), (23), (24) and (25) self-consistently, for 
the specified or desired values of the properties of sym- 
metric nuclear matter at saturation point. Further it is 
also required that the EOS so obtained has a reasonable 
nuclear incompressibility which is defined as the curva- 
ture of the energy curve at the saturation point and is 
given as. 



K = 9p', 



2 dHs/pB) 



(26) 



Incompressibility is a poorly known quantity experi- 
mentally, for the fact that some sort of theoretical mod- 
eling comes in these calculations. Apart from that, the 
other quantity with large uncertainty is the nucleon ef- 
fective mass. The wide range of values determined from 
experiments of these two quantities motivates us to an- 
alyze and study the EOS with these variations. Further 
we also need to look into issues related to some indispens- 
able elements, such as the ct— meson mass and the pion 
decay constant '/tt', while we want to achieve a proper 
framework for studying nuclear matter aspects. We recall 
that in the present work, our aim is to describe and cor- 
relate these physical quantities in a coherent and unified 
approach. 



The spontaneous breaking of chiral symmetry lends 
mass to the Hadrons and relates them to the vacuum 
expectation value (VEV) of the scalar field (xq), which 
is what is shown in eqn. (2). Immediately, what fol- 
lows from the third term in eqn. (2) is that, the VEV of 
the scalar field which has a minimum potential at /jr is 
related to the vector coupling constant C^i through the 
relation a;o = /u- = "m-uj/gLj = 1 / V Ccj ■ Thus the vec- 
tor coupling constant is explicitly constrained from the 
vacuum value of the pion decay constant. Similarly the 
scalar meson mass can be given by rria- = m VC^j/ VCa- 
The model is then sternly constrained and exude the re- 
lationship between various quantities with the VEV of 
the scalar field. 

In the present calculation, we take the value of the sat- 
uration density to be po = 0.153/to~^ [1], which agrees 
with the observed charge and mass distribution of finite 
nuclei. Saturation density implies that the pressure of the 
system is zero and the system will remain in this state 
if left undisturbed. The binding energy per nucleon is 
fixed at an empirical value B/A — m = —16.3 MeV 0]. 
With the uncertainty in the nucleon effective mass at pq, 
we calculate the parameters of the present model in the 
range Y = m* /m = (0.75 — 0.90). In order to assure the 
existence of a lower bound for the energy, we demand that 
the coefficient 'c' in the quartic scalar field term, remains 
positive. The corresponding related quantities, such as 
the pion decay constant, the scalar meson mass and the 
nuclear incompressibility are also calculated. The asym- 
metry energy coefficient is fixed at J ~ 32 MeV. The 
obtained parameters are enlisted in Table I, where the re- 
lationship between the vector coupling constant and the 
pion decay constant can be easily visualized. Stronger the 
vector coupling (repulsion), lower is the value of the chi- 
ral condensate and vice-versa. From the tabulated data, 
we find that the calculated sigma meson mass is predicted 
within (340 - 700) MeV for m*/m = (0.75 - 0.90). Al- 
though the values obtained from the analysis of neutron 
scattering off lead nuclei [1, [1^ is consistent with the 
range m* /m = (0.80 — 0.90), a lower nucleon effective 
mass is is known to reproduce the finite nuclei proper- 
ties, such as the spin-orbit effects splitting correctly [Sq . 
Also, we find that as we move to higher effective mass 
region, the incompressibility of the matter starts to fall 
or the EOS gets softer. However, within the incompress- 
ibility range of K = (200 - 300) MeV, the present model 
predicts higher nucleon effective mass. 

Nuclear matter saturation is a consequence of the in- 
terplay between the attractive (scalar) and the repulsive 
(vector) forces and hence the variation in the coupling 
strength effects other related properties as well. Fig. 
[T{A) reflects the same, where we have plotted the nuclear 
incompressibility for the evaluated parameter sets of the 
present model as a function of the nucleon effective mass. 
For better correlation between them, the corresponding 
ratio of the scalar and vector coupling is also indicated. 
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TABLE I: Parameter sets of the effective chiral model that satisfies the nuclear matter saturation properties such as binding 
energy per nucleon B /A — m — —16.3 AfeV, nucleon effective mass Y = m* /m = (0.75 — 0.90) and the asymmetry energy 
coefficient is J « 32 MeV at saturation density po = 0.153/m~^. The nucleon, the vector meson and the isovector vector meson 
masses are taken to be 939 MeV, 783 MeV and 770 MeV respectively and Co- = [ga/m^Y , = [gu/ni^Y fi^nd = {dpl'^p?' 
are the corresponding coupling constants. B = hivr? and C = c/m^ are the higher order constants in the scalar field. Other 
derived quantities such as the scalar meson mass 'jtIct', the pion decay constant '/,r' and the nuclear matter incompressibility 
(K) at Po are also given. 
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ifm^) ifm^) 


B 


C 


m„ 


Y 




K 








Um") 


(MeV) 




(MeV) 


{MeV) 


1 


5.916 


3.207 


5.060 


1.411 


1.328 


691.379 


0.75 


110.185 


1098 


2 


6.047 


3.126 


5.087 


0.822 


0.022 


675.166 


0.76 


111.601 


916 


3 


6.086 


3.031 


5.107 


0.485 


0.174 


662.642 


0.77 


113.346 


809 


4 


6.005 


2.933 


5.131 


0.582 


2.650 


656.183 


0.78 


115.238 


737 


5 


6.172 


2.825 


5.155 


-0.261 


0.606 


635.287 


0.79 


117.403 


638 


6 


6.223 


2.709 


5.178 


-0.711 


0.748 


619.585 


0.80 


119.890 


560 


7 


6.325 


2.585 


5.200 


-1.381 


0.089 


600.270 


0.81 


122.740 


491 


8 


6.405 


2.451 


5.222 


-1.990 


0.030 


580.876 


0.82 


126.039 


440 


9 


6.474 


2.323 


5.242 


-2.533 


0.300 


562.500 


0.83 


129.465 


391 


10 


6.598 


2.159 


5.265 


-3.340 


0.445 


536.838 


0.84 


134.378 


344 


11 


6.772 


1.995 


5.285 


-4.274 


0.292 


509.644 


0.85 


139.710 


303 


12 


7.022 


1.823 


5.305 


-5.414 


0.039 


478.498 


0.86 


146.131 


265 


13 


7.325 


1.642 


5.324 


-6.586 


0.571 


444.614 


0.87 


153.984 


231 


14 


7.865 


1.451 


5.343 


-8.315 


0.502 


403.303 


0.88 


163.824 


199 


15 


8.792 


1.249 


5.362 


-10.766 


0.354 


353.960 


0.89 


176.552 


168 


16 


7.942 


1.041 


5.388 


-6.908 


15.197 


339.910 


0.90 


193.437 


163 



On comparison with the incompressibility bound inferred 
from heavy ion coUision experiment (HIC)[33, we find 
that the EOS with lower nucleon effective mass is ruled 
out. The present model favors EOS for which the nucleon 
effective mass m* /ra > 0.82, i.e., the mass of the nucleon 
in the nuclear medium drops to less than w 20% of its 
mass at po- Equivalently, the agreement with the exper- 
imental flow data in the density range 2 < p/po < 4.6 
seem to favor repulsion (higher effective mass) in mat- 
ter at high density. From the plot, it can be seen that 
the EOS becomes much softer with increasing ratio of 

Ca I Cuj ■ 

Figure [TJB) shows the variation of incompressibility 
as a function of scalar meson mass obtained for vari- 
ous parameter sets. Recent experimental estimate for 
scalar meson mass m„ — 513 ± 32 MeV [l3| is compared 
with the present calculation. From the figure, we find 
that the EOS with Y = (0.84 - 0.86) (Set 10, 11 & 12; 
Table I) seems to agree with the combined constraint 
from the HIC flow data and the experimental meson mass 
range. Further, it is worth noticing that, lower the value 
of nia lower is the value of incompressibility for matter 
and vice-versa. The heavy ion collision estimate seems 
to agree with cr— meson mass within (350 — 550) MeV. 
We know that the nuclear incompressibility and the a- 
meson mass both are poorly determined quantities and 
therefore, some sort of correlation between the two will 
help to minimize the uncertainties around them. 

Figure [21 A) shows the obtained sigma mass as a func- 
tion of the vacuum value of the pion decay constant. It 



can be seen that both these physical quantities are in- 
versely proportional to each other. A lower value of f.^ 
leads to a higher value of m^. However, the experimen- 
tal bound of the pion decay constant seems to agree with 
slightly higher value of which agree with the upper 
bound of the experimental bound on rria [l^ . Precisely, 
the constraint of /jr seems to agree with ra^r ~ (560 ± 22) 
MeV. Figure [21JB) shows the nucleon effective mass Y = 
m* /m as a function of the pion decay constant. The con- 
straint of agree with EOS with Y = (0.82 - 0.84) (Set 
8, 9 & 10; Table I). However the corresponding incom- 
pressibility falls in the range K sa (344—440) MeV, which 
is on the higher side of presently acceptable bounds (Q 
- [13 )• With combined constraints such as those on nu- 
clear incompressibility inferred from HIC data [s^l, the 
limits on sigma meson mass [l3 |. the pion decay con- 
stant [l3| and the nucleon effective mass [IHl , the best fit 
from the parameters can be extracted. Thus we choose 
parameter set 9, 11 and 13 of Table I to study the cor- 
responding EOS and compare it with the experimental 
data [s^] as well as other successful relativistic mean field 
models such as NL3 parameterization pTj and the non- 
relativistic DBHF [33 calculations. 



V. EQUATION OF STATE AT T = 

The selected parameters for further study is high- 
lighted in bold fonts in Table I. The resulting energy per 
nucleon for symmetric nuclear matter is calculated for 
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FIG. 1: (Color online)(A)- Nuclear matter incompressibility 
as a function of the nucleon effective mass for the parameters 
of the present model. Also plotted is the corresponding ratio 
of scalar to vector coupling on the opposite x-axis. (B)- In- 
compressibility as a function of obtained sigma meson mass 
for the parameter sets enlisted in Table I. The upper and 
the lower limit for incompressibility inferred from Heavy Ion 
CoUision data ^ K ^ {167 - 380) MeV is shown with hori- 
zontal red lines. Recent experimental scalar meson mass limit 
{rritT = 513 ± 32) MeV [ij is depicted with blue vertical lines 



these parameters and is plotted in Fig. [3jA). For compar- 
ison, we plot the same with NL3 parameterization from 
Relativistic Mean Field (RMF) calculations [llj and also 
the non-relativistic realistic DBHF (Bonn-A) parameter- 
ization [1^. In the inset, the region of saturation density 
is magnified, where we find nice agreement within rela- 
tivistic mean- field models near pQ. Although in case of 



FIG. 2: (Color online) (A) - The scalar meson mass '771^' as 
a function of the vacuum value of the scalar field '/,r'. We 
take fn ~ 130 ± 5 MeV [l^, which is shown with red vertical 
lines. The blue lines denotes the range of the mass of ma 
from experiment [31 • (B) - Ratio of nucleon effective mass 
to bare mass is plotted as a function of 



NL3 parameter set, nuclear matter saturates at slightly 
lower density (po = 0.148 fm^^) than what we have 
taken in present calculation, the Binding energy per nu- 
cleon almost remains same {~ —16.3 MeV). In case of 
DBHF, nuclear matter saturates at still higher density. 
It is worth noticing that the incompressibility of the pa- 
rameter chosen in the present model spans within (230 
- 390) MeV, yet the resulting EOS seems to be soft at 
higher densities in comparison to that predicted by the 
NL3 parameter set which has K = 271.6 MeV. Incom- 
pressibility of nuclear matter is the measure of the degree 



of softness/ stiffness of the EOS. Conventionally, EOS 
with K < 300 MeV are eonsidered to be soft. But in the 
present case the EOS predicted by the effective model 
is relatively much softer than the NL3 parameterization 
although the value for the former is high enough. How- 
ever this can be understood, if we look at the EOS in the 
vicinity of saturation density (Inset plot). The curve of 
NL3 seems to compare well with the EOS with K — 231 
MeV below saturation density, but the energy predicted 
is much larger at higher densities. In contrast to that, 
the EOS predicted by the effective model gets softer at 
higher densities. It should be interesting to study the 
consequences of such behavior in the astrophysical con- 
text, especially on the global properties and structure of 
neutron stars at physically interesting densities (2-5 po) 

In Fig. El^B), the nucleon effective mass in the nu- 
clear medium is plotted as a function of baryon density 
up to 6/r?o- This medium mass modification of nucleon 
in nuclear medium is a consequence of the Dirac field 
and forms an essential element for the success of the rel- 
ativistic phenomenology. From the plot, it is interesting 
to see that the nucleon experiences repulsive forces in nu- 
clear matter at higher densities (ps > 2po), as a result 
of which the nucleon effective mass increases again for 
the three cases that we study presently. A careful look 
into Table I reveals the relationship between the cou- 
plings (both scalar and vector) and the resulting nucleon 
effective mass. The model predicts a higher nucleon ef- 
fective mass if the ratio of scalar to vector coupling is 
larger but the increase in m* is slower thereafter, which 
reflects the dominance of attractive force at high densi- 
ties. At saturation density, the present model results in 
much higher nucleon effective mass in comparison to the 
NL3 {m*/m = 0.60) and DBHF {m*/m = 0.678). 

Fig. mjA) displays the pressure as a function of baryon 
density up to nearly 6po for the selected parameters of 
the model for symmetric nuclear matter. The shaded re- 
gion corresponds to the experimental HIC data [33 for 
symmetric nuclear matter (SNM). Among the three the- 
oretical calculations shown, the EOS with Y ~ 0.85 & 
0.87 agree very well with the collision data. Precisely, 
the third set (K = 231 MeV) completely agree with the 
flow data in the entire density span of 2 < ps/po < 4.6. 
We now proceed to calculate the EOS of Pure Neutron 
matter (PNM) by taking the spin degeneracy 7 = 2 in 
eqn. (8) and (9). The inclusion of p— meson doesn't 
seem to affect the EOS substantially and so we refrained 
from that. In Fig. (HB), the case of pure neutron mat- 
ter (PNM) is compared with the experimental flow data. 
The experimental flow data is categorized in terms of stiff 
or soft based on whether the density dependence of the 
symmetry energy term is strong or week [4l] |. The EOS 
predicted by the present model seems to rather lie on the 
softer regime. However, the EOS with Y = 0.87, K = 231 
MeV though satisfy the combined constraint rather well, 
is not consistent with the vacuum value of the pion decay 
constant. 
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FIG. 3; (Color online) (A) - Binding energy per nucleon of 
symmetric nuclear matter plotted as a function of baryon den- 
sity up to nearly 5po- For comparison, we also plot the same 
for NL3 parameter set from the Relativistic Mean Field the- 
ory [llj as well as EOS from DBHF [s^ in the non-relativistic 
domain. The inset plot displays the curve in the vicinity of 
nuclear saturation (B) - Variation of Nucleon effective mass 
in medium as a function of total baryon density for symmet- 
ric nuclear matter of the selected parameters of the present 
model. 



VI. SUMMARY AND CONCLUSIONS 

The effective chiral model provides a natural frame- 
work to interlink the standard state properties of nuclear 
matter with the inherent fundamental constants such as 
the pion decay constant and the ct— meson mass within a 
unified approach. With this motivation, the parameters 
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FIG. 4: (Color online) Comparison of the Heavy Ion Collision 
estimate [sj with the theoretical prediction of the effective 
model, (A) for Symmetric nuclear matter (SNM) case and (B) 
for Pure Neutron matter (PNM) case. 



model is that the mass of the vector meson (m^) is gen- 
erated dynamically, as result of which the effective mass 
of the nucleon acquires a density dependence on both the 
scalar and vector fields. The interplay between this scalar 
and vector forces results in the increase of m* at « 3po 
(Fig. 3(B)), as a consequence of which the resulting EOS 
is much softer at higher densities. 

We also discussed the implication of imposing fun- 
damental constraint on the evaluated model parame- 
ters. Among various derived quantities, we find that 
the pion decay constant is experimentally well known 
quantity in comparison to the nuclear incompressibility 
and (7— meson mass, that can put stringent constraint 
on the model parameters. Employing this constraint 
(/tt = 130 ± 5 MeV) rather leave us with few options 
among the wide range of parameters enlisted in Table I, 
while the observed range of a— meson mass do not rule 
out any. Experimentally determined effective mass from 
scattering of neutron over Pb nuclei [3^ seems to go well 
with the present model. Both of them favor higher value 
for nucleon effective mass. In the present calculation we 
find that a higher nucleon effective mass is endowed with 
reasonable incompressibility too. However, the parame- 
ter that agree well (m* = 0.83to; Set 9) with the limits of 
/tt has incompressibility in the upper bound of the value 
inferred from the flow data. On a comparative analysis 
of the resulting EOS with that of the HIC data for sym- 
metric nuclear matter as well as pure neutron matter, 
parameter set with Y = 0.85; K w 300 McV seems to be 
the ideal parameterization of the present model. The re- 
sulting scalar meson mass w 510MeV, is also consis- 
tent with the experimentally observed masses |4^. [43j . 
Further, a higher value of incompressibility K « 300 
MeV is known to predict correctly the isoscalar giant 
resonance energies in medium and heavy nuclei in the 
relativistic framework [i^l . On account of the aforesaid 
arguments and constraints, the model seems to work very 
well within present approach. However, the predictabil- 
ity of the model needs to be tested at finite temperature 
and high densities. Work is in progress in this direction 
[ist - In this regard, it will also be interesting to study 
the medium effects on the underlying couplings [i^ as 
well as the on the meson masses [431 ^^'^ the pion decay 
constant [4811. 



of the model are evaluated in the mean-field ansatz by 
fixing the nuclear matter saturation properties defined 
at po and varying the nucleon effective mass in the range 
Y = m* /m = (0.75 — 0.90). Thus the resulting equa- 
tion of state not only satisfies the saturation properties 
reasonably well but also relates the various aforesaid fun- 
damental quantities with that of the vacuum value of the 
scalar field constant. One of the unique features of the 
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